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Bose-Einstein condensates of 7 Li have been limited in number due to attractive interatomic inter- 
actions. Beyond this number, the condensate undergoes collective collapse. We study theoretically 
the effect of driving low-lying collective modes of the condensate by a weak asymmetric sinusoidally 
time-dependent field. We find that driving the radial breathing mode further destabilizes the con- 
densate, while excitation of the quadrupolar surface mode causes the condensate to become more 
stable by imparting quasi-angular momentum to it. We show that a significantly larger number 
of atoms may occupy the condensate, which can then be sustained almost indefinitely. All effects 
are predicted to be clearly visible in experiments and efforts are under way for their experimental 
realization. 



Bose-Einstein condensation (BEC) in cold and dilute atomic gases 0], provide a new domain for studying nonlinear 
many-body quantum systems Consider the condensate wavefunction in the mean-field picture: First, the zero- 
point kinetic energy of the atoms tends to increase the size of the condensate wavefunction. Second, this outward 
pressure is balanced by the effect of the confining harmonic trap potential V(r). Finally, the nonlinearity enters 
through the interatomic interactions scaling as aNo\ip\ 2 , where a is the s-wave scattering length for the gas and Nq 
is the number of atoms participating in the condensate. The condensate is then well-characterized by the product 
wavefunction <f> = y / A T o-0 where the Gross-Pitaevskii (GP) equation M 
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governs the single-particle wave-function -0 of this weakly-interacting Bose condensate at zero temperature. 

Among the systems in which BEC has been observed, 7 Li is unique in that it has a negative value for a; i.e., the 
effective interaction between atoms is attractive. Such attractive interactions were initially argued to prevent BEC 
H , but it is now established that a static metastable condensate exists as long as the number of participating atoms 
iVo is less than a maximum N s Q. This is understood || to unfold as follows: As the 7 Li Bose gas is cooled below the 
critical temperature for BEC, the number of atoms in the condensate increases steadily with a kinetically determined 
'fill rate'. As Nq increases the nonlinear attractive interactions grow and the condensate becomes more localized. 
At Nq — Ag, the condensate undergoes collective collapse and there is a rapid increase in condensate density. The 
rates for inelastic collision processes, including two-body dipolar collisions and three-body recombination, are density 
dependent and increase rapidly during collapse. The atoms participating in these processes acquire large kinetic 
energies and are ejected from the condensate. Although the quantitative details of this collapse are not yet fully 
understood, experimental evidence for this process has been obtained || . 

In this Letter, we analyze theoretically the effect of weakly perturbing the harmonic trapping potential for an 
attractive condensate. We have considered general potentials V (r) = \m ufxf where r = (xi,X2,x 3 ) = (x,y, z) 
and the parameters a;, characterize the trapping potential along the three axes. For the case presented in detail 
below, lu 2 = oJq [1 + cti cos(a;t)] where u> is the frequency of the forcing and the on are its amplitudes. This corresponds 
to a spherically symmetric trap driven asymmetrically by a weak, sinusoidally time-dependent field. We employ 
a Gaussian time-dependent variational principle approximation (5| (GVA) analysis of the GP equation to solve for 
the condensate dynamics. The Gaussian ansatz is motivated by the shape of the ground-state for non-interacting 
particles in a harmonic trap and constrains the condensate wavefunction to be of the form tp{r,t) = ip x ipyip z where 

tpi = [27ro»(t)] 4 cxp j— x 2 ( 4a ^( t ) + it>i(t))^ . This ansatz when used in Eq. (|l|) yields that the condensate dynamics 
correspond to a classical effective Hamiltonian. With the transformation jl2| a,j(i) = p 2 (t) and bi(t) = ^^(t) this 
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Hamiltonian has the transparent form H var — £^ [S- -j_ _|_ m " 2 * p - ] + Zl_2^_J_ p The parameters and 11^ are 

given by (rrf) = and + PiXi) = 2/^11;. Note that (xi) = = (p^) for this wavehmction. The centroid of the 
condensate thus sits at the trap minimum and all the dynamics is restricted to changes in the condensate width pi 
in the three directions. That is, the dynamics of a quasi-particle with canonically conjugate position and momentum 
variables pi and IT characterizes the condensate. A further scaling by the length of the trap £q = \/h/ (mojo), the 
time-scale u)q , and energy Htoo yields the dimensionless Hamiltonian 



H ™ = 7 E [? + i + (! + a * cos (^))ft 2 ] - TTT~T ! 2 ! 



where /3 = 3 ^ 27r ■ The kinetics of the condensate, including the fill rate and dissipative losses, affect the dynamics 



through the time-dependence of Nq. We set the fill rate Go to a constant J7|, which is a good approximation to the 
results from the Boltzmann equation || over the time-scales considered. Further, dipolar relaxation scales as </> 2 , 
while molecular recombination is a three-body process scaling as </> 3 . In particular, with the Gaussian approximation 
for the wavefunction, we can write 

N = C ^ f —L + N ^ \ n) 

where G\ and G2 are the appropriate rate constants |8| . The Hamiltonian Eq. (||) and the kinetics Eq. ([}]) define the 
condensate dynamics completely. 

To understand these dynamics consider the fixed points given by dH var /dUi = = dH var /dpi. This yields a 
solution with = and non-zero pi corresponding to the width of a metastable condensate. This solution pp = \pp\ 
is strongly affected by the attractive hole at the origin (p = or infinite condensate density) due to the term P/(pxp2p3) 
in the Hamiltonian. The condensate is only stable when the quasi-particle avoids this attractive hole. As Nq and 
correspondingly /3 increase, the metastable minimum moves closer to the origin, becomes shallower until it becomes an 
inflection point, and finally vanishes. Correspondingly, the width pp decreases abruptly when the minimum vanishes, 
leading to collapse via the kinetics. The critical values p rn and f3 m for collapse are determined by further imposing 
the inflection condition d = 0. This gives p m = and f3 m — |5~ 5 / 4 , corresponding to iV s w 1400 atoms for 

the conditions of the Rice experiment Q . 

The normal modes of the quasi-particle dynamics around this minimum (when it exists) are the lowest-lying even 
collective excitations of the condensate |]-Jll|] . These normal mode frequencies are obtained by computing the matrices 
V,T with elements = %^Qp T . %j — oiloti an< ^ ^ en sorvm S f° r the roots of the equation det(V — fl 2 T) = 0. 



One solution is f^i = y 5 — p^ which is the frequency at which the quasi-particle vibrates radially, corresponding 

to the condensate's 'breathing' mode. The other two solutions are the degenerate solutions ^2,3 = + 2p^ 4 . 
corresponding to vibrational motion perpendicular to the radial direction. This can be visualized as being along the 
surface of a sphere of fixed radius - these are then quadrupolar surface modes |p|-|IT| of the condensate. Since pp 
decreases as a function of Nq, fli decreases and f2 2j 3 increases with increasing 7V . Also note that ^1,2, 3 = ^0 = 2 for 
pp = 1, i.e. for (3 = (non-interacting limit) in these natural units. 

We now turn to the numerical solution of this system of equations. As expected, these lowest-lying even collective 
modes of the condensate are excited by the driving. Driving the lowest such mode, a radial breathing mode, further 
destabilizes the condensate and decreases the maximum number of atoms below the static maximum N s . Remarkably, 
driving the next-highest excitations, which are quadrupolar surface modes, causes the condensate to become more 
stable. That is, for certain frequencies, atoms occupy the condensate in numbers significantly greater than TVs and 
the condensate is sustained for correspondingly longer times. We have studied in particular the frequency-dependent 
value of the maximum number of atoms N u that can be sustained in the presence of the driving. We find that this 
frequency response is broad and qualitatively robust, with similar features for wide range of tested parameters and 
many alternative configurations for both the confining and perturbing fields. An example of these responses is shown in 
Fig. (1), along with a theoretical estimate iV an g as explained below. In obtaining it, we use = 0.02 = —2a! = — 2a 2 
which represents the effect of a Helmholtz field oriented along the axis of an Ioffe-Pritchard magnetic trap. The 
response is shown as the ratio N m = N u /N s . The broad stabilization effect is clear in this figure; further, it can 
be seen that the curve has some interesting structure, including in particular a signature dip corresponding to a 
parametric resonance between the radial and quadrupolar modes as explained below. The stabilization effect is so 
strong that in some cases it is no longer obvious whether and when a collapse occurs. It is possible for Nq to grow so 
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large that the loss rate equals the fill rate, even for relatively large p. The values in Fig. (1) are therefore calculated 
using a 'worst case' definition of the collapse, as the Nq value at which the inelastic decay terms in Eq. (|J) exceed 
the growth rate Go- Although it is unclear whether the physics of the true dynamical situation in these cases is fully 
captured by the simple kinetics of Eq. (|j|), the stabilization is nonetheless remarkable. 

To understand the detailed structure of Fig. (1), consider the case of driving the condensate with u> < CIq. Note 
that the dynamics of the quasi-particle are always restricted to the p\ = P2 plane by our choice of driving. At an 
arbitrary value of w, the driving is not initially resonant with cither of the excitations and there is a negligible response 
from the condensate. As Nq increases with time, however, f2i(/3) will ultimately equal any us < 51o. The breathing 
mode is excited as a result and the quasi-particle oscillates along the radial direction. As this oscillation increases 
in amplitude, the quasi-particle is driven into the attractive hole and the condensate collapses. The net effect is to 
decrease N m (u>), with the impact being greatest for us slightly less than Slo- The more interesting case is for us > Qq 
where again, the driving is not initially resonant. As No increases, however, the quadrupolar surface modes come 
into resonance. The quasi-particle then oscillates with increasingly greater amplitude in a direction perpendicular 
to the radial direction, still in the p\ = P2 plane. An angular momentum vector lq(t) — p x II pointing in the q 
direction can be associated with this motion, where q = -^(1,-1,0) is the normal to the p\ = p 2 plane. This 
angular momentum vector reverses orientation |l3] during the oscillation, with its magnitude going to zero at the 
turning points of the oscillation. However, the mean-square value of this angular momentum increases as the surface 
mode acquires increasing energy. By virtue of the energy in this oscillation, the quasi-particle avoids the chasm of 
the attractive 'hole', and the condensate is stabilized, avoiding collapse even when the attractive interactions are 
significant. The above implies that the stabilizing effect may be quantitatively estimated by adding to the kinetic 
energy pressure terms an angular momentum term ^ a where I 2 is the time-average of Z? (t) . This yields an estimated 

increase in the maximum number of atoms as A ang = ^4j=S = (1 + ^-) 5 ^ 4 ', this estimate is plotted in Fig. (1) with 

I 2 computed directly from the dynamics. It is clear from comparing the curves that this simple angular momentum 
argument captures the essential features of the stabilization. In particular, the region of discrepancy is precisely where 
the naive kinetics and collapse criterion render suspect the value of A w shown, as argued above. 

This resonant driving is easier to sustain for frequencies close to J7o since f22,3 has a weaker dependence on Nq 
in that regime. As iVo increases, the condensate ultimately falls out of resonance with the driving. Consider for 
example the condensate radius p{t) — \p{t)\: It is initially unaffected by the driving and decreases slowly as in the 
static case; this is followed by a resonantly driven increase that ultimately saturates. An example of this mechanism 
is shown in Fig. (2a) for us = 2.3 (corresponding to a driving frequency of s» 308-ffz for the Rice experiment). In 
Fig. (2b) we show a numerical solution for the GP equation corresponding to the situation in Fig. (2a); the two 
curves agree qualitatively. We note here that the numerical solution to the nonlinear Schrodinger equation in the 
GP form for cylindrical symmetry is extremely difficult and a full scan of frequencies is computationally prohibitive. 
However, simulations including filling and inelastic decay similar to those in the literature [fl3| have been performed 
to qualtitatively validate the GVA results. 

As us increases, the resonance occurs at increasingly later times and for decreasing windows of time, and hence the 
greatest stabilization happens for us slightly greater than Qq. There is no such stabilization for us > us 2 = >/l2 = 3.46; 
this critical frequency is obtained by substituting the maximum value for the equilibrium p m into the expression for 
^2,3 _ i.e., the driving cannot resonate with the condensate modes if the metastable minimum does not exist in the first 
place. This cut-off can be seen clearly in Fig. (1). Finally, we point out the significant dip in the curves at u) = u>i = y8. 
This frequency corresponds to the situation where fl2,s(0) = 20i(/3) i.e. for pj 4 = 3. In this case, the driving and 
the surface modes parametrically excite the radial mode through the 1 : 2 resonance between the frequencies. Thus, 
the condensate is destabilized even for a frequency greater than Qq . Similar resonances account for the other detailed 
structure in the frequency response curves. All these features should be clearly visible experimentally, since the 
experimental resolution is of the order of 50 condensate atoms || , and the predicted features exceed this resolution 
significantly. Moreover, the details of the curves are sensitive to the precise models and parameters used. Comparison 
of experimental results with these theoretical predictions will therefore help a) define the regime of validity of a given 
model and b) improve our understanding of the details of the rich interplay between kinetics and non-linear many- 
body quantum dynamics in the attractive condensate. Efforts are under way for experimental realization of these 
phenomena. 

The instability and finite lifetimes of attractive condensates have been a constraint in the push to understand 
these systems. It is expected that the significantly increased stability through weak driving will open the way to 
further novel experimental phenomena and interesting theory and ultimately a deeper understanding of this nonlinear 
quantum regime. 
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FIG. 1. Stabilization of an attractive condensate by weak driving as measured by the change in the maximum number of 
atoms participating in the condensate as a function of the driving frequency. The ratio iV m between the frequency-dependent 
maximum number of atoms N u and the maximum value iV s for the static condensate is shown. The figure is truncated at 
iV m ?s 4.8 due to the excessive (t > 100s) amount of time taken to reach this value. An estimated value A r ang for N m obtained 
from an angular momentum analysis is also shown. See text for details. 

FIG. 2. (a) Time-dependence of the radially-averaged width of the condensate p(t). There are rapid oscillations, as expected, 
on the time-scale of the driving (308-ffz) that cannot be easily resolved on the scale of this figure, (b) As in (a) except as 
computed with a numerical solution of the GP equation with added dissipation. See text for details. 
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